Low density ferromagnetism in Extended Hubbard ( t-t'-U ) Model 

P. A. Sreeram, Suresh G. Mishra 
Institute of Physics, Bhubaneswar 751005, India 

We study the existence of ferromagnetism in one dimensional Hubbard Model with 
on site interaction and nearest and next nearest neighbor hopping. Using the Hubbard 
I approximation, a self consistent equation for < n a >, the particle density with spin 
a is obtained. We find that ferromagnetism exists over a wide range of values of the 
next nearest neighbor hopping element t' when on-site interaction U is large. The phase 
diagram as a function of t' /t and the electron density < n > has been obtained by 
numerically solving the self consistent equation. It is found that the maximum density 
at which ferromagnetism can occur, peaks around t'/t ~ —0.53. We show that this 
phenomenon is a result of the peculiar way in which the t! term modifies the density of 
states. 

PACS numbers :05.60.+W,72.10.Bg,67.57.Hi 

For last thirty years various techniques have been used in the study of the existence of ferromagnetism 
within the Hubbard Model. Within the t-matrix approximation Kanamori Q had shown, that, for a given 
on-site repulsion U there exists a critical density p c below which no ferromagnetism can exist. Mean Field 
theory, however, allows the existence of ferromagnetism as long as the Stoner Criterion Up(ef) > 1 is 
satisfied, where /?(e/) is the density of states at Fermi energy, ej. In his first paper (henceforth, we shall 
refer to this as Hubbard-I) Hubbard himself had asserted that, if the density of states is peaked 
away from the center of the band, then the ferromagnetic phase can exist for low density of electrons. 
The Lieb-Mattis theorem ||, however, precludes any possibility of existence of ferromagnetism in the 
1-dimensional Hubbard Model for any arbitrary filling. 

On the other hand, in certain limiting cases, e.g. the Nagaoka point |||| , at U = oo and a one hole 
in an otherwise half filled band it is possible to show the existence of ferromagnetism. Another example 
where saturated ferromagnetism has been shown to exist in Hubbard Model is the case of Flat Band 
Ferromagnetism ||. In this mechanism, the ground state of the Hamiltonian H = Ho + Hi„ t is at the 
same time the ground state of Ho and Hj n t separately. This model can be mapped on to a degenerate 
band Hubbard Model. 

Recently Miiller-Hartman has shown that it is possible to get ferromagnetism in 1-D Hubbard 
Model at low electron densities, if one introduces next-nearest neighbor hopping (t 1 ), besides the nearest 
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neighbor hopping matrix element t, for the condition t' /t < —0.25. This work has been extended by 
Pieri et. al. J8| using the Gutzwiller approach. In the present work, we study the t — t' — U model using 
the Hubbard I approximation. We find Miiller-Hartman result reminiscent of the comment made in the 
Hubbard first paper. The introduction of the t' term brings more weightage to the lower edge of the 
density of states, moreover it also takes care of the limitation concerning the Lieb Mattis theorem. We 
now give some detail of our work. 

The model under consideration is given by the Hamiltonian, 

H = -tJ2( c l c i+i° + ct+^Ci) - t'J2(clc i+2 * + c] +2a c ia ) + uY d n i ^n ii (l) 

ia i i 

where C\ a (Ci a ) is the creation (annihilation) operator of an electron with spin a at the site i, ni a = 
C\ a Ci a and t,t' and U have the usual meaning, as defined earlier. The noninteracting Hamiltonian can 
be diagonalized in the momentum space giving, 

H (k) = J2^ClC ka (2) 

kcr 

where takes the form, 

e k = -2tcos(k) -2t'cos(2k), (3) 
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FIG. 1. Free electron dispersion curve for t /t = 0.0,-0.25 and -0.5. 

for one dimension. The dispersion is parabolic at low k and develops a double well structure as the value 
of t' jt reduced below - 0.25 (Fig. 1). It is shown below that the existence of ferromagnetism in the system 
is closely related to this particular structure of the dispersion. 
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The Greens Function, G a (k,oj) —<< Cy, C], >>, for the Hamiltonian given in Eq. 1, can be calculated 
using the Memory Function Technique H; it is given by, 



G a (k,u) = [u> -e k - U < n_ a > (1 + — — F ) + irf\ 1 

LU — U 



(4) 



upto the first order in the memory kernel, where, U = U(l— < n- a >), and is given by Eq. 3. The 
number density < n a > is related to the Greens function by, 



1 

27r J-oo u 



(5) 



k 

We thus have a self consistent equation for < n a > . The phase diagrams obtained by solving Eq. 5 are 
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FIG. 2. Phase diagram in the U-< n > plane for t' /% = 0.0, -0.1, -0.25, -0.52 and -0.8. The ferromagnetic phase 
exists below the curve corresponding to each value of t'/t. 



shown in Figs. 2 and 3. In Fig. 2 we show the phase diagram as a function of electron density < n > 
and the on-site repulsion U, for various values of the next nearest neighbor hoping parameter t'. From 
the figure, it is obvious that for very low value of U, the region of ferromagnetic phase is very small. As 
U increases, the critical density of electrons < n > c up to which ferromagnetism can exist increases, but 
this critical density saturates beyond a certain value of U. Thus, for any given value of t'/t, there exists a 
maximum value of < n > c beyond which ferromagnetism cannot exist, even if U is increased indefinitely. 
The dependence on t'/t is also clear from this figure. For large U ( when < n > c has saturated ), there 
exists a optimum value of t'/t (« -0.52) which gives maximum < n > c . This becomes very clear from the 
phase diagram shown in Fig-3. However, at low values of U, this optimum value of t'/t shifts towards « 
-0.25. 

Fig. 3 shows the phase diagram in the < n >-t' plane for U—10. The value of U has been such that 
the critical value of < n > c has reached the saturated value for large U. Thus this phase diagram is a 
reflection of the large U behaviour. We have chosen this regime since the Hubbard-I solution tends to 
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overemphasize correlation and hence is quite unreliable at low values of U. The phase diagram shows that 
the density of electrons below which ferromagnetism can occur increases as t'/t is reduced below zero, 
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FIG. 3. Phase diagram in the t' /t-< n > plane for (7=10. The shaded region denotes the ferromagnetic regime. 
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FIG. 4. Density of states for 17=10 and t'/t = (a) 0.0, (b) -0.1, (c) -0.25, (d) -0.52, (c) -0.8. The figure shows 
only the Lower Hubbard Band since it is the only band which is relevant to our discussion. The free electron 
density of states looks very similar. 



passes through a maximum at t'/t w -0.52 and then starts decreasing. In order to understand the phase 
diagrams shown in Figs 2 and 3 we have studied the behaviour of the density of states as a function of 
t'/t. The explanation for ferromagnetism in the t-t'-U Model lies in the fact that the lower edge peak in 
the density of states is far away from the center of the band. This gives rise to ferromagnetism as argued 
in Hubbard-I. As t'/t is reduced from zero, initially up to t'/t =-0.25, this peak grows in height (Fig. 



4), thus giving rise to larger and larger densities at which ferromagnetism can occur. As t' jt reduces 
towards -0.52, the height of the lower peak starts reducing again but it begins to develop a spread, which 
can accommodate larger density of electrons. Beyond t' jt rj -0.52, this spread begins to move towards 
the center of the band and hence, the critical density < n > c is reduced. Thus, at t' ft sa -0.52, the area 
within the lower peak has an optimum value, thus giving a maximum value of < n > c . 

In conclusion, the critical density < n > c below which ferromagnetism is possible in the t-t'-U Model 
is governed by the density of states argument given in Hubbard-I. There are, however, some issues, 
which must be discussed here. As mentioned earlier, the Lieb-Mattis Theorem states that there can 
be no ferromagnetism for t' = ( which corresponds to the usual Hubbard Model) in one dimension. 
However, we find that within the Hubbard-I approximation, ferromagnetism occurs at t' jt = 0, at very 
low densities. This is essentially because of the approximation involved, whereby, although correlation is 
taken into account beyond the Hatree Fock approximation, the solution, in essence, remains a mean field 
solution. However, in spite of the mean field nature of the phase diagram at t' jt ^ gives an indication 
of the trend of the actual phase diagram which can be obtained by more accurate and exact calculations, 
work on which is in progress. 
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